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Ab-initio calculations for the structural properties of Zr — Nb alloys at different values of the 
niobium concentration are done at zero temperature. Cases of BCC and HCP structures are con- 
sidered. Optimal values of the lattice constant together with structural relation c/a for the HCP 
crystals are obtained. Critical value for the niobium concentration which corresponds to the struc- 
tural transformation HCP — >■ BCC at zero temperature in alloys Zr — Nb is defined. Densities of 
states for two different structures with niobium concentrations 12.5% and 25% which have HCP and 
BCC structures, accordingly, are studied. 

PACS numbers: 61.50.Lt, 61.72.Bb, 61.72.jj, 61.72.S-, 71.15.Dx, 71.15.Nc 



I. INTRODUCTION 

Investigation of Zirconium alloy properties is an actual 
task during last 60-80 years in Nuclear Physics and Power 
Engineering and has very special practical significance. 
It is due to Zirconium and Zirconium-Niobium alloys are 
widely use in nuclear power industry as materials for the 
fuel elements, fuel assemblies and other constructions of 
nuclear reactors [iHll- It is known, that Zirconium has 
two configurations: a-phase, which is characterize by the 
Hexagon Close Package (HCP) lattice and stable at not 
too high temperature (< 863 '^C), and /3-phase with Body 
Centered Cubic (BCC) lattice at high temperature. Nio- 
bium has BCC lattice. Zirconium-Niobium alloys which 
characterize by a-phase are used as construction material 
of the active zone of nuclear reactors because of success- 
ful a successful mix of nuclear, mechanical and corrosive 
properties 

Alloys Zr — xNb with a high percentage of Niobium 
{x > 5%) have a high capacity for hydrogenation, re- 
sulting in deterioration of corrosion properties. That is 
why for the reactor metallurgy Zirconium- Niobium alloys 
with small Niobium concentration (< 5%) are useful. It 
is because in the nuclear reactor Hydrogen diffuses in- 
side an alloy due to corrosion processes. When the Hy- 
drogen concentration becomes more than solubility limit, 
it stands out as a fragile hydride phase, that limits the 
usage of an alloy in the reactor. Hence, the study of 
Zr — Nb alloys with small and high Niobium concen- 
tration, together with processes of hydride formation in 
Zirconium-Niobium alloys is an actual task. Investiga- 
tion of the hydride phase formation processes in alloy 
Zr — Nb that contains 20% and 40% of Niobium was 
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made in dll]. Tu study such process in alloys with high 
Niobium concentration the plane strain theory was used 
13. In % diffusion coefficients of Zr and Nb atoms in 
BCC alloys with Niobium concentration 5.5%, 16.3% and 
28.1% were studied. It was found, that Zirconium and 
Niobium atoms diffuse in alloy according to the mono- 
vacancy mechanism. The last is defined by the specifica- 
tion of lattice dynamics. Diffusion of Hafnium and Nio- 
bium atoms in coarse-grained BCC alloy Zr — Nb with 
19% Niobium was studied in Q- It was defined, that 
diffusion coefficient of Niobium is less the Hafnium one. 
Corrosion behavior of alloy Zr — Nb with Niobium con- 
centration variation from 0.02% to 20% was considered 
in Q . It was shown, that corrosion rate becomes smaller 
with decrease in both bulk of the /3-phase and Niobium 
concentration i a-phase. Authors defined that the most 
important factor that reduces corrosion rate is the pre- 
cipitate of the /3 — Nb particles and the corresponding 
decreasing of Niobium in the Zirconium matrix. Experi- 
mental investigations of /3-phase formation with Niobium 
concentration 20-^88% at irradiation were made in [S-J^] . 
Experimental results for the structural properties a- and 
/3-phases in Zr — Nb alloys with varying in Niobium con- 
centration at high temperature (~ 700K) are presented 
in [T^ - [l6j . Structural properties of Zr — 12.5%Nb were 
theoretically considered in [13]. However, the question 
about the critical value of Niobium concentration when 
the transition from a-phase toward /3-phase (structural 
transformation HCP — > BCC) in Zirconium-Niobium al- 
loys is steel open. 

In this work we are aim to study structural properties 
of Zirconium- Niobium alloys with different Niobium con- 
centration at low (zero) temperature from the first prin- 
ciples. Model structures of Zr — Nb alloys with BCC and 
HCP lattices will be considered. The main task of this ar- 
ticle is to define optimal lattice constants and in the case 
of HCP structures, optimal structural relation c/a for dif- 
ferent Niobium concentration in the alloy. We shall de- 
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fine the critical value of Niobium concentration that cor- 
responds to the structural transformation HCP — BCC 
at low (zero) temperature. The quantitative comparison 
with known experimental data will be done. 

The work is organized in the following manner. In Sec- 
tion 2 we present models for studied structures of Zr—Nb 
alloys and research methods description. In Section 3 
we discuss structural properties and obtain optimal lat- 
tice parameters for studied structures. Analysis of total 
energy and density of states is done in Section 4. We 
conclude in the last Section. 



II. MODELS OF THE STUDYING 
STRUCTURES AND RESEARCH METHODS 

It is known, that Zirconium at low temperatures (< 
863°C) is characterized by /3-phase (HCP structure), 
whereas Niobium has BCC structure. Hence it is nat- 
urally to expect, that alloys Zr — xNb will have HCP 
structure at low Niobium concentration x whereas BCC 
structure will be stable for large x. As far as Zirconium 
and Niobium are different only for one electron and are 
characterized by close package structures one can expect 
that Niobium atoms in the HCP Zirconium lattice, or 
Zirconium atoms in the BCC Niobium lattice are substi- 
tutional atoms. In Ref.[18] authors have studied super- 
cell of 48 atoms of HCP Zirconium which has a vacancy 
and Niobium atom in interstitial position placed near 
the vacancy. It was shown, that during relaxation pro- 
cesses Niobium atom takes a vacancy position. The last 
confirms that atoms of Zr or Nb in Zr — Nb alloys are 
substitutional atoms. 

In this article we will consider alloys Zr — xNb with 
BCC and HCP lattices, where the Niobium concentra- 
tion X takes values x = 3.125%, 6.25%, 12.5%, 25.0%, 
50.0% and x = 75.0% for both types of lattices, and 
additionally x = 87.5%, 93.75% and x = 96.875% for 
BCC alloys. In such a case, structures with Niobium 
concentration 3.125% and 96.875% are characterized by 
the largest number of atoms in the unit cell, 32. Hence, 
for the given values of the Niobium concentration x in the 
alloy Zr — xNb with the unit cell of 32 atoms, a number 
of Niobium atoms can be : 1, 2, 4, 8, 16, 24, 28, 30 and 
31, respectively. 

Model structures of BCC alloys with 2, 16 and 30 Nio- 
bium atoms in the unit cell of 32 atoms are characterized 
by the 221 symmetry space group; others structures have 
the 123 group of symmetry. All being studied structures 
of Zr — Nb alloys with HCP lattice are characterized 
by the 187 space group. The number of nonequivalent 
atoms of Zirconium and Niobium, that gives a number of 
symmetry operations for each model of BCC and HCP 
structure is shown in Table [Hi 

All calculations of structural and electronic properties 
for the model alloys, shown in table HIl were done in the 
framework of density functional theory (DFT) [l^ by 
using linearized augmented plane wave (LAPW) method. 
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TABLE I: The number of Niobium atoms in the unit cell 
of 32 atoms {if^Nbuc) and numbers of nonequivalent atoms 
of Zirconium {ij^NEAzr) and Niobium (i^NEANb) for each 
model structure with BCC and HCP lattices. 

which is implemented in software packages Wien2k j20l |. 
This method self-consistently includes basic and valence 
electrons and is widely used for the calculations of the 
band structure of solids [20| . 

The Muffin-tin radii {Rmt) of both Zr and Nb were 
taken to be 2.1 atomic units for all calculations. The basis 
function was expand up to Rmt x kmax — 7, where Rmt 
is the smallest radius of the muffin-tin (MT) spheres and 
kmax is the maximal value of the reciprocal lattice vec- 
tors. Wave function expansion inside the atomic spheres 
was done until Imax = 10. Integration over the Brillouin 
zone was performed by using 1000 fc-points. The criterion 
to stop the iteration process was converge for the total 
energy and charge of the crystal to less than 0.0001 Ryd 
and 0.001 e~^, respectively. All calculations were done 
using generalized gradient approximation (GGA) with 
parametrization (PBE) [21,] . 

III. STRUCTURAL RELATIONS 

In this section the main results for the structural prop- 
erties of studied structures will be shown. To define op- 
timal values of the lattice constants for each structure of 
Zr — Nb alloy and structural relation c/a for the HCP 
crystals the optimization procedure was used. 

A. BCC lattice 

For alloys Zr — xNb with BCC lattice the optimization 
procedure was made in the standard manner. For each 
value of Niobium concentration x we have defined the 
total energy for the unit cell of the crystal with varying 
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in the unit cell volume. Obtained data were approxi- 
mated by the equation of state [121 ■ The minimum on 
the dependence of the total energy of unit cell versus unit 
cell volume corresponds to the optimal value of the unit 
cell volume. Following the standard definition for the 
volume of the unit cell of the BCC-crystals Vbcc = afcc 
optimal value for the lattice constant abcc was defined. 
The dependence of the optimal lattice constant abcc on 
the Niobium concentration x for BCC alloys Zr — xNb 
is shown in FigHJ 




at. %Nb 

FIG. 1: Dependence of the lattice parameter abcc for the BCC 
alloys Zr ~ Nb versus Niobium concentration. 

Here, obtained theoretical results are shown with the 
help of filled squares. Empty and filled trian gles cor- 
respond to the experimental data (see Refs. \23\ and 
[14| . respectively). The known value for the lattice con- 
stant of pure BCC Niobium is shown by symbol x . From 
Fig ID it is seen that with an increase in the Niobium con- 
centration in the BCC alloy Zr — xNb the lattice con- 
stant decreases according to the linear law. This result 
is consistent with the well known approximate empiri- 
cal rule, called in metallurgy Vegard's law [23|. Accord- 
ing to this rule at constant temperature a linear depen- 
dence between the crystal lattice parameter of an alloy 
and the concentrations of the constituent elements is re- 
alized. Thus for the studied BCC alloy NbxZri^x, where 
Zirconium atoms are substitutional ones in the BCC Nio- 
bium lattice, the relation between lattice parameters for 
the pure Zirconium and Niobium and its alloy is given 
as follows: of'zrNb = ^ ' + (1 ^ ^) ■ ■ Here value 
for the lattice constant of pure BCC Niobium is known, 
a^g = 3.3004 A, whereas a|=;= = 3.5634 A— fitting pa- 
rameter (at low temperatures Zirconium is characterized 
by the HCP lattice). 

It should be noted that in the interval of Niobium con- 
centration X from 10% to 50% in the BCC alloy Zr — xNb 
experimental data are slightly exceed obtained theoret- 
ical ones. It is natural, due to the theoretical calcula- 
tions were performed at low (zero) temperature, while 
the presented exp erimental results obtained at evaluated 
temperatures [ij, [1^ . It was experimentally shown, that 
in the crystal Zr — xNb at (~ TOOK) BCC phase is real- 
ized at evaluated and large values for the Niobium con- 



centration, i.e. at.%Nb > 10% [H,!!!]. In FiglUwe also 
show values for lattice constant of BCC alloy at small 
Niobium concentrations. In the next Section making an 
analysis of the total energy values of the crystals at low 
(zero) temperature we shall define the minimal value of 
the Niobium concentration in the alloy Zr — xNb, when 
crystal will characterize by BCC lattice. 



B. HCP lattice 

Now let us consider alloys Zr — xNb with HCP lattice. 
The optimization procedure allowing to find optimal lat- 
tice constants ahcp and Chcp for such crystals is more dif- 
ficult then in the previous case for BCC lattice. It is due 
to that for the HCP alloy one should define a global mini- 
mum of the total energy of a crystal as a function of both 
the volume of unit cell and the structural relation c/a. 
To this end we have proceeded in the following manner. 
For each value of Niobium concentration a test value of 
structural relation c/a was fixed. The standard optimiza- 
tion procedure, like in the case of BCC lattice was used 
to define the optimal volume of unit cell corresponding 
to the minimum of the total energy as a function of the 
unit cell volume, El^-^{{c/a)i), where i = 1...N, N is 
a number of test values for the structural relation c/a. 
Obtained empiric data as dependence Emmie/ a) were 
approximated by a functional dependence. Minimum of 
this dependence corresponds to the optimal value of the 
structural relation {c/a)opt for each studied value of Nio- 
bium concentration x in the HCP alloy Zr — xNb. Next, 
for defined {c/a)"^^^ the standard optimization proce- 
dure was made to define a global minimum of the total 
energy and the optimal unit cell volume, respectively. 
Using the standard definition for the volume of the HCP 
unit cell Vhcp = (c/a)^^^ a^^^ sin(27r/3) the optimal lat- 
tice constant a^cp was defined for each studied value of 
the Niobium concentration in the alloy. Obtained results 
for the lattice constant aucp and the structural relation 
c/a as functions on the Niobium concentration x in the 
HCP alloy Zr — xNb are shown in FigsH^ and ^jp, re- 
spectively. 

In Figures [5^ and [2Jd obtained theoretical results are 
shown by circles. By symbol x values for the pure HCP 
Zirconium are shown. First, let us consider dependence of 
the lattice constant ahcp versus Niobium concentration, 
shown in Figl2^. It is seen, that increase in the Niobium 
concentration x in the HCP alloy Zr — xNb leads to de- 
crease in the lattice constant ahcp, like in the case of the 
BCC crystals. However, in contrast to the previous case, 
the falling down dependence ahcp{at.%Nb) differs from 
the linear for the shown interval of the Niobium concen- 
tration values. At low values of Niobium concentrations 
{at.%Nb < 20%) calculated data for the lattice constant 
a-hcp lie on the linear dependence, but for large at.%Nb 
one has deviation from the linear law. Thereby, empiric 
Vegard's rule for the studied HCP alloys can be applied 
in the case of low values of Niobium concentration x in al- 
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FIG. 2: Dependencies of the lattice parameter ahcp (a) and 
structural relation c/a (b) for the HOP alloys Zr — Nb versus 
Niobium concentration. 



loys Zr — xNb. In the insertion in Figl2^ we have shown 
obtained theoretical results for ahcp at low at.%Nb to- 
gether with experimental data presented by triangles and 
squares from and [1^ , respectively. One can see good 
quantitatively correspondence with experimental results 
for small values of Niobium concentration. 

Next, let us consider dependence of the structural re- 
lation c/a as a function on the Niobium concentration 
X for the HCP alloys Zr — xNb shown in Figl^b- Here, 
marking is made in the same way, as in Fig|2^ for the lat- 
tice constant. It is seen, that an increase in the Niobium 
concentration in the alloy leads to a grow of the struc- 
tural relation c/a. It means, that the unit cell becomes 
more prolongate in z-direction. As for the lattice con- 
stant ahcp for the structural relation c/a the falling down 
dependence versus Niobium concentration differs from 
the linear law; the last can be applied only for small val- 
ues of at.%Nb. Comparison of results obtained for HCP 
Zr — Nb alloy with experimental data from Refs.fU, [2^ 
for structural relation c/a at at.%Nb < 10% gives small 
deviation, which is less than 1% (see insertion in Figl^t). 

Experimental investigations at T ~ 7Q0K indicate that 
HCP phase in Zr — Nb alloys is realized only if at.%Nb < 
6.5% [3, HI]. Next, analyzing values of the total energy 
the maximal Niobium concentration value in alloy Zr — 
Nb, when alloy still characterizes by HCP lattice at low 
(zero) temperature will be defined. 



IV. TOTAL ENERGY AND DENSITY OF 
STATES 

Previously we have calculated optimal values for lat- 
tice constant for alloys Zr — xNb with BCC and HCP 
lattices at different values of Niobium concentration x. 
Due to at normal conditions pure Zirconium character- 
ized by the HCP lattice but pure Niobium has BCC lat- 
tice one can expect that at fixed temperature there is a 
critical value for Niobium concentration Xc correspond- 
ing to the structural transformation HCP -> BCC. It 
means, that for alloy Zr — xNb ai x < Xc energetically 
more favorable will be HCP structure, whereas at x > Xc 
BCC structure will be realized. Experimental investiga- 
tions of Zirconium-Niobium alloys at high temperatures 
(T ~ 700K) show that in such conditions the critical 
Niobium concentration lies in the interval from 7 to 10 % 

To define the critical value of the Niobium concentra- 
tion Xc at low (zero) temperature we shall analyze the 
total energies of all studied crystals with both HCP and 
BCC lattices and determine the energetically most favor- 
able structure (HCP or BCC) for each model. To that 
end we proceed in the following manner. For fixed Nio- 
bium concentration x using obtained optimal lattice con- 
stant afjcc, a-hcp and, in the case of HCP lattice, structure 
relation c/a we make a full cycle of ab-initio calculations 
to define the total energy for the unit cell. The unit cell 
contains 32 atoms for each structure due to the maximal 
number of atoms (32) in the unit cell have structures 
Zr — 3.125%Nb (1 atom of Niobium and 31 atoms of Zir- 
conium) and Zr — 96.875%iV6 (31 atoms of Niobium and 
1 atom of Zirconium). Next, we define the difference be- 
tween total energies for the unit cell of 32 atoms of HCP 
and BCC lattice, as: 

AEUat.%Nb) = Er:,{at.%Nb)\^^^-EZ{at.%Nb)\,^^. 

Hence, the value AEtot{at.%Nb) defines the type of en- 
ergetically most favorable lattice for the alloy Zr — xNb 
at fixed value of Niobium concentration x as follows: (i) 



if El^,iat.%Nb)\ 



hep 



< 



E^^\(at.%Nb)\ 



and as a result 



AEtotiat.%Nb) < 0, then alloy Zr - xNb has HCP lat- 
tice; (ii) at Er,%{at.%Nb\^^ > E^S,{at.%Nb)\,^^, that 
gives AEtot{at.%Nb) > 0, the BCC lattice is realized. 
Therefore, the critical value of the Niobium concentra- 
tion Xc (or {at.%Nb)c) that defines the structural trans- 
formation HCP — ?> BCC at fixed temperature (low/zero 
temperature in the studied case) can be defined from 
the condition of energetic equivalence of two lattices, i.e. 
AEtot{at.%Nb) = 0. Obtained results for the difference 
of the total energies of the unit cell of 32 atoms versus 
Niobium concentration in alloy Zr — Nb are shown in 
Figm 

It is seen that AEtot changes a sign from minus to 
plus in the vicinity of the critical value of Niobium con- 
centration {at.%Nb)c — 20%. Therefore, at low (zero) 
temperature the alloy Zr — Nb with the Niobium concen- 
tration less than 20% is characterize by the HCP lattice 
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FIG. 3: A difference in the total energy of HCP and BCC 
structures AEtot{cit-%Nb) in eV for the unit cell of 32 atoms 
as a function of the Niobium concentration in percents for 
alloys Zr - Nb. 



{AEtot < 0); values for the difference AEtot are shown 
by circles. In the case of large Niobium concentration, 
more than 20%, the Zirconium-Niobium alloy has BCC 
lattice (squares in the FiglS]). Obtained result gives an 
explanation about the deviation from the linear depen- 
dence of the obtained results for the lattice constant Uhcp 
of the HCP alloys versus Niobium concentration for large 
at.%Nb (see FigUi). 

As was pointed out before, experimental investiga- 
tions for the Zirconium-Niobium alloys at high temper- 
atures 700K) show, that the BCC phase is realized 
if at.%Nb > 10%, whereas the HCP phase can be ob- 
served ifat.%Nb < 6.5% [Hill,!!!]. Thus, theoretically 
obtained result from ab-initio calculation for the criti- 
cal value of Niobium concentration in the Zr — Nb al- 
loy {at.%Nb)c — 20%, that corresponds to the structural 
transformation HCP — > BCC is about twice of the exper- 
imental one. An increase of the critical value of the Nio- 
bium concentration with a decrease in the temperature 
seems to be natural. With the temperature increase one 
gets that thermal fluctuations become large in the vicin- 
ity of the lattice knots. It leads to the destruction of the 
more complex HCP structure at smaller Niobium concen- 
tration. Hence one can expect a falling down dependence 
of the critical Niobium concentration {at.%Nb)c in alloy 
Zr — Nb versus temperature. 

As far as obtained critical Niobium concentration at 
low (zero) temperature is {at.%Nb)c — 20%, next we 
perform calculations of the density of states (DOS) for 
two structures: Zr - 12.5%Nb and Zr - 25%Nb, which 
are characterize by different energetically more favorable 
lattices, HCP and BCC, respectively. Obtained results 
are shown in Fig|4l Here, solid curves correspond to al- 
loys with HCP lattices, whereas dash ones relate to al- 
loys with BCC lattices. It is seen, that for the alloy 
Zr — 12.5%A^6 the amplitude of the main peaks in the 
density of states is larger for the BCC lattice (see FigH^), 
whereas for the structure Zr — 25%Nb situation is quite 
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FIG. 4: Dependencies for densities of states (DOS) versus 
energy measured from the Fermi energy for: (a) alloy Zr — 
12.5%Nb; (b) alloy Zr - 25%Nb in the case of HCP lattice 
(solid curves) and BCC lattice (dash curves). 



different: larger DOS corresponds to the alloy with HCP 
lattice (see FigHjD). These results confirm that at low 
(zero) temperature alloy Zr — 12.5%Nb is characterize 
by the HCP lattice, whereas alloy Zr - 25%Nb has BCC 
lattice. 



V. CONCLUSIONS 

We have studied structural properties of the 
Zirconium-Niobium alloys with different Niobium con- 
centrations in the framework of ab-initio calculation, us- 
ing the density functional theory (DFT), linearized aug- 
mented plane wave (LAPW) method and generalized gra- 
dient approximation (GGA) with parametrization (PBE) 
at low (zero) temperature. Niobium concentration was 
varied from 3.125% to 96.875%. Models with BCC and 
HCP lattices for each considered Niobium concentration 
in alloy were studied. Using the optimization procedure 
to find a global minimum of the total energy of a crys- 
tal the optimal values for the lattice constants, and in 
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the case of HCP lattice, structural relation c/ a were ob- 
tained. 

We have found, that introduction of Zirconium atoms, 
as substitutional ones into the BCC crystal of Niobium 
leads to increase in the lattice constant. Obtained de- 
pendence of lattice constant versus Zirconium concentra- 
tion in the BCC Niobium alloy is in good correspondence 
with empiric Vegard's law and is quantitatively consis- 
tent with known experimental data. Introduction of sub- 
stitutional Niobium atoms into HCP Zirconium crystal 



results to the decrease in the lattice constant and 
result to the increase in the structural relation c/a. The 
correspondence with the Vegard law, here, is realized only 
at small Niobium concentration (less then 20%), whereas 
obtained results are consistent with experimental ones. 
We have defined that at low (zero) temperature the crit- 
ical value of the Niobium concentration in alloy Zr — Nb, 
which defines the structural transformation HCP BCC 
is {at.%Nb)c 20%. 
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